Quantum Phase Transition of Ground-state Entanglement in a Heisenberg Spin Chain 

Simulated in an NMR Quantum Computer 
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Using an NMR quantum computer, we experimentally simulate the quantum phase transition of 
a Heisenberg spin chain. The Hamiltonian is generated by a multiple pulse sequence, the nuclear 
spin system is prepared in its (pseudo-pure) ground state and the effective Hamiltonian varied in 
such a way that the Heisenberg chain is taken from a product state to an entangled state and finally 
to a different product state. 
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Quantum mechanical systems are known to undergo 
phase transitions at zero temperature when a suitable 
control parameters in its Hamiltonian is varied At 
the critical point where the quantum phase transition 
(QPT) occurs, the ground state of the system undergoes 
a qualitative change in some of its properties Oster- 
loh et al. [2| showed that in a class of one-dimensional 
magnetic systems, the QPT is associated with a change of 
entanglement, and that the entanglement shows scaling 
behavior in the vicinity of the transition point. This be- 
havior was discussed in detail for the Heisenberg model 
3] and for the Hubbard model p|. It is believed that 
the ground-state entanglement also plays a crucial role in 
other QPTs, like the change of conductivity in the Mott 
insulator-superfluid transition |5| and the quantum Hall 
effect Many of the relevant features, like the transi- 
tion from a simple product state to a strongly entangled 
state, occur over a wide range of parameters and persist 
for infinite systems as well as for systems with as few as 
two spins 00- These systems, especially the Heisen- 
berg spin model, are central both to condensed-matter 
physics and to quantum information theory. In quantum 
information processing, the Heisenberg exchange interac- 
tion has been shown to provide a universal set of gates 
0, [ljj an d m quantum communication, information can 
be propagated through a Heisenberg spin chain . 

While some Heisenberg models can be solved analyti- 
cally, others can only be simulated numerically. Like for 
other quantum systems, such simulations are extremely 
inefficient if the system contains more than 10-20 spins. 
It was therefore suggested that such simulations could be 
more efficiently performed on a quantum computer [l2T |. 
In this Letter, we discuss the simulation of a Heisenberg 
spin chain by an NMR quantum computer. By vary- 
ing the strength of the magnetic field, we take the sys- 
tem, which is in the quantum mechanical ground state, 
through the QPT and measure the change in entangle- 
ment by quantum state tomography. The NMR tech- 
niques that we use here are closely related to earlier 



work where they were used to demonstrate quantum algo- 
rithms, quantum error correction, quantum simulation, 
quantum teleportation and more(see, e.g., Ref. [13j and 
references cited therein). 

The simplest system that exhibits this behavior con- 
sists of two spins coupled by the Ising interaction 



Jiaial, J i > 0, 



(1) 



where a\ are the Pauli operators, ui z a magnetic field 
strength, and Jj is a spin-spin coupling constant. 

In the range — Jj < tu z < Ji, the ground state of this 
system is two-fold degenerate. To avoid this complica- 
tion, we add a small transverse magnetic field. The re- 
sulting Hamiltonian is 



f(*i+o2) + (2) 



which is nondegenerate. The transverse field will always 
be kept small, \u x \ <C \u> z \, \ Ji\- 

A symmetry-adapted basis that is an eigenba- 
sis for vanishing transverse field (uj x = 0) is 

{in),i* + },iu},i*-}}, with |*±) = -l(iti) ± in)) 

and ||) and ||) the spin-down (m = — |) and spin up 
(to = +i) states. Furthermore, it is convenient to define 
dimensionless field strengths g x = and g z = |j. 

Since the ground state of this system is always one 
of the triplet states, and transitions to the singlet state 
are symmetry-forbidden, we can reduce our system of 
interest to the triplet states. For small transverse fields, 
g x <C 1, the longitudinal field g z determines the ground 
state 




9z < -1 

-Kg z <1 
9z > 1 



(3) 



g z = ±1 are therefore quantum critical points, where the 
ground state changes from the ferromagnetically ordered 



2 



high field states to the entangled, antiferromagnetic low- 
field state. 

For the full system, including the transverse field, the 
eigenstates and eigenvalues of the three-state system are 



\^ = ^ t+m+ ^-^ ITT) 

+ ^^l* + ) + IU}),(^l,2,3) 



and 



£1,2 = J/6,2 = J/(l - 2r cos(# T tt/3))/3 
£3 = = Jj(2r cos 0+l)/3, 

where Mi are normalization constants, r 
V3(5^+5, 2 ) + l, and = | arccos( 4(18g '-"'- 2) ). 



(4) 



(5) 




FIG. 1: (a) Energy level diagram for the two-spin Heisenberg 
Ising model for g x — (dashed lines) and g x = 0.129 (solid 
lines), (b) Probability amplitudes of the ground state for 
9x = 0.129. 



where we assumed that the rf field strength is the same 
for both qubits. 

The target Hamiltonian (0 can be created as an av- 
erage Hamiltonian by concatenating small flip angle rf 
pulses with short periods of free evolution, e~ tHr ~ 
fT Pe -iH NM Rr prS a^ wn ere t p is the pulse duration and 
the length of the free evolution period. The rcsult- 



,-tJ?, 



' prec 

ing effective Hamiltonian matches the target Hamilto- 

- Ml, 

'prec — j 12 



nian if uj\ = lo 2 l = 
While this approximation is correct to first order in r, 
the symmetrized sequence 

(Hrft _ (HNMR,Tprec) ~ {H r f, y ) (8) 

generates the desired evolution to second order in r. Fig- 
ure^Jshows the sequence of rf pulses required to generate 
this evolution. 
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FIG. 2: Sequence of rf pulses applied to both spins to simu- 
late the target Hamiltonian. The boxes represent pulses that 
induce rotations around the x-axes of the rotating frame, the 
separations between them the free precession periods. The 
index m for the different periods runs from to M. 



Figure ^ shows numerical values for the energies and 
the coefficients of the ground state as a function of the 
longitudinal field strength g z = jjj. The right hand 
side shows clearly that at strong fields (\g z \ > 1), the 
ground state is a product state, while it corresponds to 
the entangled state \^ + ) for weak fields. 

To observe the system undergoing the QPT, we simu- 
late it on an NMR quantum computer, where the quan- 
tum spins a 1 are represented by nuclear spins and the 
Hamiltonian of the Heisenberg chain is mapped into 
an effective Hamiltonian generated by a sequence of radio 
frequency (rf) pulses acting on the nuclear spin system. 

The natural Hamiltonian of our two-qubit system is 



Ji 



(6) 



The 



1,2 



represent the Larmor frequencies of the two 



qubits and Jyi the spin-spin coupling constant. In addi- 
tion to this static Hamiltonian, we use rf pulses to drive 
the dynamics of the system. In the usual rotating co- 
ordinate system, the effect of rf pulses can be written 
as 



H 



rf 



-FT\Px +< J x) 



(7) 



To prepare the system in the grou nd state, we use the 
technique of pseudo-pure states |1J|: we prepare a den- 



sity operator p pp (ipi) 



trJT) 



1 + a|V'i)(V'i|- Here, 1 is 



the unity operator and a a small constant of the order of 
10~ 5 . To measure the order parameter (entanglement), 
we apply quantum state tomography 1 1 5| . The system 
can then be taken through the phase transition by adi- 
abatically changing the magnetic field g z of the effective 
Hamiltonian, which acts as a control parameter. 

To ensure that the system always stays in the ground 
state, the variation of the control parameter has to be 

QMt)|^(t)) 



sufficiently slow, so that the condition 



= (*)"£!(*) 



< 1 



is fulfilled, where the index e refers to the excited states 
Choosing g z as the control parameter, we write the 
adiabaticity condition as 



dg z 
dt 



« Jf X - Jj 



(&(*)- &(t)) 2 



(Mt)\W-\Mt)) 



(9) 



where the dimensionless parameter \ quantifies the sen- 
sitivity to the control parameter g z and we have concen- 
trated on the first excited state IV'2}, which is the critical 
one for transitions from the ground state. 
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FIG. 3: Adiabatic magnetic field sweeps g z (t). The solid line 



was calculated for constant adiabaticity parameter 



dg z I dt 



(s 



equ. ©) for a transverse field g x — 0.129; the circles represent 
the values obtained from the numerical optimisation of the 
discretised scan. 



Equation @ defines the optimal sweep of the control 

<'</= 



parameter g z (t), with the scan speed ^§f- oc \- Figure 
13 shows the required time dependence of the magnetic 
field. The resulting transfer is therefore highest for a 
given scan time or the scan time minimised for a required 
adiabaticity. 

The experimental implementation generates an effec- 
tive Hamiltonian that is constant for a time r (see equ. 
©). For this stepwise approximation, the duration of 
each time step has to be chosen such that (i) the time 
is short enough that the average Hamiltonian approxi- 
mation holds and (ii) the adiabaticity criterion remains 
valid. While this calls for many short steps, there is also a 
lower limit for the duration of each step, which is dictated 
by experimental aspects: switching transients, which are 
not taken into account in the Hamiltonian of equ. ijHJ, 
tend to generate errors that increase with the number of 
cycles. 

We used a numerical optimisation procedure to deter- 
mine the optimal sequence of Hamiltonians, taking the 
full level structure into account. Choosing a hyperbolic 
sine as the functional form, we optimised its parameters 
and found the optimised discrete scan represented by the 
circles in Fig. |31 

To determine the optimal number of steps, we used the 
same numerical simulation, keeping the functional depen- 
dence g z vs. t constant, but increasing the number of 
steps. The results are summarized in Fig. 0] which plots 
the lowest fidelity encountered during each scan against 
the number of steps taken in the simulation. The fidelity 
is calculated as the overlap of the state with the ground 
state at the relevant position. The simulation shows also 
the effect of decoherence, which reduces the achievable 
fidelity if the total duration of the scan becomes compa- 
rable to the decoherence time. The model that we used 
to take the effect of decoherence into account is similar 
to that of Vandersypen et al . 

For the experimental implementation, we used the 



FIG. 4: Numerical simulation of the minimum fidelities dur- 
ing the adiabatic passage vs. number of steps with (*) and 
without (O) decoherence effects. 



13 C and X H spins of 13 C_labeled chloroform (both spins 
1/2). The relatively large spin-spin coupling constant of 
4 12 =214.94Hz makes this molecule well suited for this ex- 
periment. The chloroform was diluted in acetone-dg and 
experiments were carried out at room temperature on a 
Bruker DRX-500 MHz spectrometer. The pseudo-pure 
initial state p PP (TT) was generated by spatial averaging 
[Tij . The fidelity of this state preparation was checked by 
quantum state tomography and found to be better than 
0.99. 

The adiabatic scan through the QPT was achieved by 
shifting the rf frequencies of both channels by the same 
amount after each period. Using the sweep g z {t) shown in 
Fig. El the offset was changed from g z = — 3 to g z = +3 
in 60 steps. The evolution of the system during the scan 
was checked by performing a complete quantum state to- 
mography after every second step during the experiment. 

As a quantitative measure of the QPT, we used the 
concurrence as the order parameter, which is related to 
" the entanglement of formation" [l|| and ranges form 
(no entanglement) to 1 (maximum entanglement). For 
this purpose, we calculated the concurrence C from the 
tomographically reconstructed deviation density matri- 
ces p as C(p) — max{\i — X2 — A3 — A4,0}, where 
Aj(i = 1,2,3,4) are the square roots of the eigenvalues of 
/ 3 ( cr y cr y)/ 3 *( cr y cr y) m decreasing order. 

Figure Efa) shows the measured concurrence C cxp = 
C(/°cxp) as individual points and compares them with 
the theoretical values Cth- Both data sets clearly show 
the expected QPTs near the critical points g z = ±1. The 
entangled ground state for \g z \ < 1 is characterised by a 
concurrence close to 1, while the high field states are only 
weakly entangled (the entanglement vanishes for g x =0). 

The experimentally determined concurrence remains 
below ~ 0.75, significantly less than the theoretical val- 
ues. To verify that this deviation is due to decoherence, 
we simulated the experiment, taking into account the de- 
tails of the pulse sequence as well as the effect of decoher- 
ence. We obtained good agreement between theoretical 
and experimental data if we assumed a total decoherence 
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FIG. 5: (Color online) (a) Measured fidelity F cxp (+) and con- 
currence Ccxp (*) compared to the concurrence calculated for 
an ideal scan Cth (solid line) and the simulated concurrence 
Cdec (v) an( i fidelity Fd ec (O) when decoherence is taken into 
account, (b) Measured values (*) of the two-spin correlaror 
(°"zO"z) compared to the theoretical (solid line) and simulated 
values with decoherence (v). 



time of 130 ms, which is slightly longer than the 110 ms 
scan time used in the experiment. Fig[5Ja) shows the 
simulated values of the concurrence as triangles; their 
evolution during the scan is quite similar to that of the 
experimental data points. 

To assess the quality of the adiabatic scan, we also 
determined the fidelities -F 0X p = F(p cxp ) from the to- 
mographically reconstructed density operators. The fi- 
delities, which are shown at the top of Fig[5Ja), deviate 
from unity when the system passes through the critical 
points and shows some overall decrease due to decoher- 
ence. Again, the simulated fidelities agree remarkably 
well with the experimental values. 

As a second order parameter, we also determined the 
two-spin correlation .1] (o^erf) = Tr(p exp ala^), which 
are shown in Fig|3fb). As expected, the system is ferro- 
magnetically ordered ((cr^crf) = +1) at high fields, but 
turns to an antiferromagnetic state ((er^cr^) = —1) at low 
fields between the two quantum critical points. Compar- 
ing Fig. |5fa) with (b), the concurrence has the similar 
behavior to the two-spin correlation. 

In conclusion, we have discussed an experimental quan- 
tum simulation of a quantum phase transition in a 
Heisenberg spin chain. Heisenberg spin chains, which 
have been investigated in detail in solid state physics, 
play an important role in a number of proposed solid 
state quantum computers. During the course of the sim- 
ulation, the system ground state changes from a classical 
product state to an entangled state and back to another 



product state. Like in many other proposed quantum 
simulations, this system had to be swept adiabatically 
through the relevant parameter space by properly vary- 
ing a Hamiltonian parameter. The techniques developed 
here may also be useful for other types of adiabatic quan- 
tum computing which have been proved to be effective 
for NP-hard problems [l^, e -g-j the classically NP-hard 
ground state search. Furthermore, the adiabatic passage 
can provide a novel method to create entanglement [2£|. 
The simulation can be extended to other types of Heisen- 
berg spin chains, e.g., Heisenberg XY or XYZ models 
etc.; work in this direction is under way. 
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